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Fisher's zeros as boundary of renormalization group flows in complex coupling spaces 
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We propose new methods to extend the renormahzation group transformation to complex couphng 
spaces. We argue that the Fisher's zeros are located at the boundary of the complex basin of 
attraction of infra-red fixed points. We support this picture with numerical calculations at finite 
volume for two-dimensional 0{N) models in the large- A*' limit and the hierarchical Ising model. We 
present numerical evidence that, as the volume increases, the Fisher's zeros of 4-dimensional pure 
gauge SU{2) lattice gauge theory with a Wilson action, stabilize at a distance larger than 0.15 from 
the real axis in the complex /? = 4,/g^ plane. We discuss the implications for proofs of confinement 
and searches for nontrivial infra-red fixed points in models beyond the standard model. 
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The renormalization group (RG) method has played 
a prominent role in advancing our understanding of 
strongly interacting and strongly correlated systems. A 
question of central importance in this context is to find 
out if the physical spectrum of gauge theories, sigma 
models or Hubbard models contains massless particles, 
such as gauge or Nambu-Goldstone bosons, or if a mass 
gap is generated dynamically. For asymptotically free 
theories, this question can be addressed by studying the 
(marginally relevant) RG fiows coming out of the weakly 
coupled fixed point. If these RG flows end at the strongly 
coupled fixed point, the correlations among local observ- 
ables decay exponentially with the separation (mass gap) 
and Wilson loops decay exponentially with their area 
(confinement). Showing rigorously that this statement 
correct is a possible strategy [l| to prove confinement in 
non-abelian gauge theories on the lattice. 

Confinement can be lost by either introducing a finite 
temperature or enough matter fields to modify the run- 
ning of the coupling constant. In particular, if enough 
species of fermions are added without spoiling asymptotic 
freedom, one may expect a nontrivial infrared (IR) fixed 
point with conformal symmetry j^. Recently, there has 
been a renewed interest in this possibility to build mod- 
els for possible new physics beyond the standard model. 
Various extensions of QCD have been studied [3,]. Other 
scenarios involving multiple confinement/ deconfinement 
transitions were also proposed in Refs. [J|. Establishing 
the existence of an IR fixed point is often controversial 
and it would be desirable to find criteria independent of 
the RG method used. 

In this Letter, we show that considerable insight on 
these questions can be obtained by extending the RG 
fiows in complex coupling space. We provide empiri- 
cal evidence that the global behavior of the complex RG 
fiows can be determined by simply calculating the com- 
plex zeros of the partition function in the complex cou- 
pling plane (Fisher's zeros) and related singular points. 
In the large volume limit, we argue that the basin of at- 
traction of the strongly coupled (confining) fixed point is 



delimited by Fisher's zeros. The complex conjugated pair 
of zeros closest to the real axis can be seen as a "gate" 
controlling the complex flows between the weakly and 
strongly coupled fixed points. Monitoring the position of 
these zeros as the volume increases can provide a way to 
decide if the theory is confining or not. For SU{2) pure 
gauge theory this is a challenging task Q especially when 
we increase the volume. In the following, we report the 
first numerically stable calculations of Fisher zeros on a 
G"' lattice. Details will be provided elsewhere Q. 

In the study of flows or differential equations, it is often 
enlightening to consider their complexification. Recent 
RG studies Q discuss the loss of conformal invariance 
and disappearance of fixed points in the complex plane 
when a parameter is varied beyond some critical value. In 
addition, analytic continuation of a theory into the com- 
plex coupling plane is an essential tool to understand the 
large-order behavior of perturbative series. It has been 
used to explain why perturbative series have a zero ra- 
dius of convergence [8| and to determine accurately the 
growth of the perturbative series^ 9, 10]. For lattice mod- 
els with compact integration [ll], [I^l , there is a change 
but not a loss of vacuum when the real positive coupling 
5q changes sign. This may explain the apparent power 
growth behavior of weak coupling series obtained with 
stochastic perturbation theory |13l - [l5| . A complex RG 
perspective on these questions would be very desirable. 

In the following, we consider three types of models 
known for their absence of phase transition: the non- 
linear 0{N) sigma model on a square lattice in the large- 
N limit, the two-dimensional Ising hierarchical model 
and SU{2) lattice gauge theory in 4 dimensions with a 
Wilson action. We use generic notations which can be 
used interchangeably for the 3 models. (3 denotes the 
inverse 't Hooft coupling l/{gQN) for the 0{N) model, 
the inverse temperature for the hierarchical model, and 
4/5q for the SU{2) gauge theory. These three models 
are discussed in more detail in Refs. [12,|lg, [l^l respec- 
tively. We use the notation a for the lattice spacing, 
me for the mass gap and M = aniQ its dimension- 



less form. /? should not be confused with the Callan- 
Symanzik /3 function which will be denoted Pes- In gen- 
eral, M^dp/dAP oc Pes /do, with a model dependent 
positive constant of proportionality. For asymptotically 
free theories, in the small a limit we have 



P{M^) ~ A + B In(l/Af2) 



(1) 



with B > 0, proportional to minus the first coefficient of 
the /Scs-function. We will also use the generic spectral 
decomposition [ITI, lla] of the partition function: 



Z = 



dSn{S)e-f^^ 



(2) 



where S is the total action (or energy). We call n{S) 
the density of states. We use the notation ln(n(5))/7V = 
/{S/Af) with TV the number of lattice sites for spin mod- 
els and the number of plaquettes for gauge models. 

Simple one-dimensional RG flows can be generated 
when the /3 dependence on the lattice spacing a at fixed 
physical mass gap mQ is known. Beyond the celebrated 
asymptotic scaling regime of Eq. ([T]), local polynomial 
parametrizations of In(M^) are known in SU{3) lattice 
gauge theory jl9| . but their analytical continuation is 
only valid in a small region. For the 0{N) models at 
infinite volume and infinite N, we have a close form ex- 
pression [2(| valid for any real positive M^: 



P{M 



-f2x 



d^k 



_J_A'27rrPHkuk2) + M\ 



(3) 



with P^(fci, ^2) = 2(2-cos(fci)-cos(fc2)). In the limit of 
small M^, we have the asymptotic form of Eq. ^ with 
B ^ l/(47r). 

The mapping can be analytically continued to the cut 
complex M^ plane with a cut running from -8 to 0. 
The image of this cut plane in the complex /? plane is 
an asymptoticall y cr oss-shaped region partially shown in 
Fig. [T] In Ref. [i3| it was argued that the Fisher's ze- 
ros should lay outside of the this crossed shaped region. 
Complex RG flows can be obtained by increasing the lat- 
tice spacing (kept real) with a fixed complex rriQ. Fig. 
[T] shows the RG flows for 11 initial values of M^ taken 
on a small circle around the origin and then multiplied 
repeatedly by a factor 2. The flows stay inside the image 
of the [—8,0] cut. By taking initial values very close to 
the real negative axis, it is possible to follow closely this 
boundary. 

Lattice simulations are performed at finite volume and 
it is important to understand the modification of the sim- 
ple flow picture provided above when the volume is finite. 
For L X L lattices with periodic boundary conditions, the 
integral in Eq. (|31) is replaced by an average over the L^ 
momenta coming in units of 27r/L. The mapping /3(M^) 
becomes a rational function and its inversion requires a 
Riemann surface with a number of sheets scaling like L^. 
The cuts (in the p plane) end at the images /3(M|j„g ) of 
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FIG. 1: Top: Infinite L RG flows (arrows). The blending 
small crosses (blue on-line) are the /3 images of two lines of 
points located very close above and below the [—8,0] cut. 
Bottom: same procedure and initial conditions but for L = 
32; the crosses are the images of the singular points. The 
image of the CSP described in the text appear as two large 
filled circles. 



values of P can be interpreted as the complex zeros of 
the non-perturbative /Jcs-function [12| . If we consider 
the image of a straight line coming out of the origin in 
the complex M^ plane and making an angle 9 with the 
positive real axis, its image may wrap around a certain 
number of the P{M^^^ ). A detailed analysis ^21)] shows 
that it only occurs when \9\ > tt/2. The complex conju- 
gate images with the smallest 6 play an essential role in 
the understanding of the global properties of the flows. 
In practice, their 9 is very close to ±7r/2 and we call them 
the closest singular points (CSP). The situation is illus- 
trated in the bottom part of Fig. [l]for L = 32. The pro- 
cedure is identical to the infinite volume case, but some 
results are strinkingly different. For very small |M^|, the 



pole at in the finite volume sum that replaces the inte- 
gral in Eq. ([3]), dominates and replaces the logarithmic 
divergence at infinite volume. Consequently, the image of 
a small circle around the origin in the M^ plane is a large 
circle in the /3 plane. There are 288 singular points. The 
flows corresponding to \9\ < 7r/2 go between the images 
of the CSP (large filled circles) from the right. The flows 
corresponding to 16*1 > 7r/2 go between the images of the 
CSP from the left, wrap around the images of the CSP 
and eventually end up at zero. In summary, the global 
properties of the RG flows are controlled by the CSP. In 
the large volume limit, the real part of the image of the 
CSP goes to infinity and the imaginary part stabilizes at 
(7r/2)(l/(47r)) = 1/8 as can be inferred from Eq. (P). 
This stabilization implies that the complex RG fiows on 
the positive real axis and the neighboring complex fiows 
can reach the strongly coupled fixed point without ob- 
struction. 

Complex RG flows can also be calculated by extend- 
ing 2-lattice matching methods to the case of complex 
/3. In the following, we use a slightly modified version of 
Ref. [22, [23 . The idea is to consider large distance ob- 
servables that can be calculated on lattices with different 
sizes in lattice spacing units but equal physical sizes. The 
large distance behavior is probed by calculating the cor- 
relations between two large neighbor blocks B and NB 
of physical size \B\. The ratio of the block volume to the 
total physical volume V is the same for both lattices. In 
order to bypass the determination of the field rescaling, 
we consider the ratio 



HM matching n=4 and n=5 



D\ — 



R{f5,V/a^) = 



iT,xeB'l^x)iT,yeBM) 



(4) 



In Ref. [23,], the whole ratio was averaged. Here, we 
use blocked observables that depend linearly on the orig- 
inal variables, the average at complex /3 can be defined 
by reweighting configurations at real l3. A discrete RG 
transformation mapping /3 into /3' while the lattice spac- 
ing changes from a to ba is obtained by requiring the 
matching R{P, V/a^) = R{j3' , V/iba)^). When ^ is com- 
plex, there are typically many /3' . In special cases, the 
matching condition reduces to polynomial equations in 
which the multivaluedness can be addressed systemati- 
cally. For practical purposes, one would like to be able 
to use Newton's method to construct the RG flows. This 
works if there is only one /?' solution close to /3. For the 
two spin models considered here, we found out that un- 
less the RG flow get near the Fisher's zeros, the distance 
1/3 — /3' I singles out one /3' unambiguously. The situation 
is illustrated for the hierarchical model in Fig. [21 We 
required R{/3,2^) — R{l3',2*) using the exactly calcula- 
ble probability distribution for blocks covering half the 
volume. 29 initial /3 were chosen on a line with constant 
Re/S = 5. The 8 trajectories passing by the Fisher's 
zeros led to ambiguous choices of /3' and are not dis- 
played. More detail on this method will be provided in 
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FIG. 2: Unambiguous RG flows for the hierarchical model in 
the complex /3 plane obtained by the two lattice method. The 
crosses and open boxes are at the Fishers zeros for 2^ and 2^ 
sites. 



Refs. (2l[,[2J]. 

We now introduce a generic method to restrict the re- 
gion where Fisher's zeros can be found. Using the den- 
sity of states, it is clear that the contributions at fixed /3 
should come from a small region near the saddle point. 



n{S)e 



-l3Afs _ „A/'(/(s)-/3s) 



(5) 



-A/'(/(so) + (l/2)/"(so)(s-so)'-H...) 



with s = S/N and /'(so) = /3. As long as Ref"{so) < 0, 
the distribution becomes Gaussian in the infinite volume 
because if we define Aff"{so){s — sq)^ = v"^ as the normal 
variable, higher order in s — sq in the exponential will be 
suppressed by negative powers of J\f. Gaussian distri- 
butions have no complex zeros [25|. and consequently, 
we could look for the level curve Ref"{so) = as the 
boundary of the re gion where Fisher's zeros may appear. 
In Fig. 1 of Ref. [2a| the regions where Ref"{so) > 
are depicted as narrow "tongues" coming vertically to- 
ward the real axis. In the U(l) case, a conjugate pair 
pinches the real axis, but for SU{2) a finite gap remains 
present. This suggests that the Fisher's zeros of these 
models should appear on approximately vertical linear 
structures. For SU{2), the imaginary part of Fisher's 
are too large to use simple reweighting methods [5]. By 
using Chebyshev interpolation for /(s) and monitoring 
the numerical stability of the integrals with the residue 
theorem [I^l , it is possible to obtain reasonably stable re- 
sults [y] that confirm this picture (see Fig. [3]) . Unlike the 
C/(l) case, the imaginary part of the lowest zeros does not 
decrease as the volume increases, but their linear density 
increases at a rate compatible with L"*. 

For 0{N) models, it is possible to write close form ex- 
pressions of the partition function in the approximation 
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FIG. 3: Images of the zeros of /"(s) in the /? plane (open 
symbols) and Fisher's zeros (filled symbols) for 4^ (squares) 
and 6* (circles) lattices. 



(justified in the large- TV limit) where we only keep the 
zero mode of the auxiliary field enforcing the constraint 
(/>.(/> — 1. For L and N not too large, it is possible to 
use the residue theorem to calculate exactly the integral. 



From this exact expression, one can calculate the Fisher's 
zero and the density of states which happens to be piece- 
wise polynomial. Using the exact form of the density of 
state, we can calculate the zeros of /". As in the case of 
the SU{2) gauge theory, approximately vertical lines of 
zeros appear above the singular points of the two map- 
pings discussed above. Their linear density increases at 
a rate compatible with L~^ |2]| . 

In summary, we have shown with examples that RG 
flows in one real coupling can be extended to the complex 
coupling plane. As the volume increases, the stabiliza- 
tion of the Fisher's zeros away from the real axis allows 
the complex flows to reach the strongly coupled fixed 
point. Two- lattice matching methods can be extended 
to the complex plane by reweighting existing gauge or 
spin configurations (as long as the imaginary part of (3 is 
not too large). We plan to apply this method to decide 
if extensions of QCD confine or not. 
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